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1 Introduction and Motivation 



Hadronic light-cone distribution amplitudes (DAs) of leading twist play an essential role in 
the QCD description of hard exclusive processes. They can be obtained from the hadron's 
Bethe-Salpeter equation by integration over the transverse momentum distribution, but 
keeping track of the longitudinal momentum fraction u: 



DAs enter the amplitudes of processes to which coUinear factorisation theorems apply and 
which notably include "classical" applications like the EM form factor of the pion or the 
77*-7r transition form factor, which were first discussed in the seminal papers by Brodsky, 
Lepage and others [1]. More recently, collinear factorisation has been shown to apply, to 
leading order in an expansion in l/rrib, also to a large class of nonleptonic B decays [2], 
which has opened a new and exciting area of applications of meson DAs: nonleptonic B 
decays, and in particular CP asymmetries in these decays, are currently being studied at 
the B factories BaBar and Belle and are expected to yield essential information about the 
pattern of CP violation and potential sources of flavour violation beyond the SM. DAs 
also enter crucially SCET, the soft-collinear effective theory [3], which aims to provide a 
unified theoretical framework for the factorisation of both hard- collinear effects, relevant 
for all hard exclusive QCD processes, and soft effects, which occur in processes involving 
heavy mesons. 

For light mesons, on which we shall concentrate in this paper, leading twist DAs can 
be interpreted as probability amplitudes of finding the meson in a state with a minimum 
number of Fock constituents and at small transverse separation which provides an ultra- 
violet cut-off. The dependence on the this cut-off is given by Brodsky-Lcpage evolution 
equations and can be calculated in perturbative QCD, while the DAs at a certain low scale 
provide the necessary nonperturbative input for a rigorous QCD treatment of exclusive 
reactions with large momentum transfer. 

As exclusive QCD processes come with much smaller cross-sections or branching ra- 
tios than inclusive processes, they are not as well studied experimentally as their inclusive 
counterparts. The overall normalisations of DAs are given by local hadronic matrix el- 
ements, essentially decay constants, which are partly accessible experimentally, partly 
have to be calculated from theory. Information on the shape of DAs comes mostly from 
theory and has been the subject of numerous studies within various nonperturbative ap- 
proaches. There are a few, mostly exploratory, studies of the second momentum of the 
pion DA in lattice QCD, Ref. [4], but the results are not yet refined enough to be relevant 
for phenomenology. Most of the existing information on light meson DAs comes from 
QCD sum rules, whose application to the study of moments of DAs has been pioneered 
by Chernyak and Zhitnitsky. Their review on the "Asymptotic Behavior Of Exclusive 
Processes In QCD" , Ref. [5] , is still a very useful source of information, despite having 
been written 20 years ago. Later analyses switched from the calculation of moments to 
that of Gegenbauer moments which, as we shall discuss in Sec. 2, are more appropriate 
for the study of DAs. The leading-twist DA of the pion was analysed in Ref. [6], while 
those of the p were studied in Ref. [7]. The first paper to study SU(3) breaking effects in 
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moments was Ref . [8] , while the first two Gegenbauer moments of the K* and $ meson 
DAs were calculated in Ref. [9]. 

The aim of the present paper is to provide a careful analysis of SU(3) breaking effects 
in leading-twist K and K* DAs, using QCD sum rules. The motivation for this study 
is the need for a more accurate assessment of these effects than presented in the papers 
quoted above. This need is driven by the fact that DAs are crucial ingredients in the 
analysis of B decays of type B nK etc., which are presently being measured at the 
B factories. In the future, accurate information on SU(3) breaking effects in DAs will 
also be required for assessing the reliability of phenomenological methods for relating, via 
U-spin symmetry, amplitudes of Bd decays to those of Bg decays which will be measured 
at the Tevatron and the LHC. SU(3) breaking enters these processes in essentially three 
different ways: 

• via an asymmetry of the DAs under the exchange of quark and antiquark; 

• via SU(3) breaking in the symmetric parts of the DA, in particular the overall 
normalisation; 

• via form factors, e.g. F^^'^ vs. F^^^ . 

Although in the present paper we will concentrate on the first two manifestations of SU(3) 
breaking, our results can be used immediately for an update of B ^ K, K* form factors 
using QCD sum rules on the light-cone (cf. [10] for the present state of the art). 

Our paper is organised as follows: Sec. 2 contains the definition of all relevant DAs as 
well as a discussion of their behaviour under a change of renormalisation scale. In Sec. 3 
we discuss QCD sum rules for the DAs, obtain numerical results and compare with the 
results of other studies. Section 4 presents a summary and our conclusions. Issues of a 
more technical nature are discussed in the appendices. 

2 General Framework 
2.1 Definitions 

We define the light-cone DAs via matrix elements of quark-antiquark gauge-invariant 
nonlocal operators at light-like separations with = Q [5] . For definiteness we consider 
the K^*^- meson distributions; the DAs of the neutral mesons containing an s quark just 
involve a trivial isospin factor l/-\/2 in the overall normalisation. For mesons containing 
an s antiquark, one has (f){sq){u) — (t>{qs){^ — u). The complete set of distributions to 
leading-twist accuracy involves three DAs (we use the notation z — z^^^ for arbitrary 
4- vectors z): 

{Q\u{z)z-i^[z,Q]smK-{q)) = ifK{q - z) [' due-'^'^'^-'^Kiu) , 

Jo 

{0\u{z)z[z,0]siO)\K*-{q,X)) = ie^^h)flmK* /' diie--^^-)^^, 

Jo 
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{0\u{z)a^4^,0]siO)\K*-{q,X)) 



with the Wilson-hne 



[z, 0] = Pexp 



da z'^An{az) 



due-'''^'i-'Ui{u), (2.1) 



inserted between quark fields to render the matrix elements gauge-invariant. In the above 
definitions, ei'' is the polarization vector of a vector meson with polarisation A; there 
are two DAs for vector mesons, corresponding to longitudinal and transverse polarisation, 
respectively. The integration variable u is the meson momentum fraction carried by the 
quark, u = 1 — u the momentum fraction carried by the antiquark. The normalisation 
constants Jk are defined by the local limit of Eqs. (2.1) and chosen in such a way that 



du 4>{u) 



(2.2) 



for all the three distributions 



(Pk, (Pk- 



2.2 Conformal Expansion and Renormalisation 

The conformal expansion of light-cone distribution amplitudes is analogous to the partial 
wave expansion of wave functions in standard quantum mechanics. In conformal expan- 
sion, the invariance of massless QCD under conformal transformations is the equivalent 
of rotational symmetry in quantum mechanics, where, for spherically symmetric poten- 
tials, the partial wave decomposition serves to separate angular degrees of freedom from 
radial ones. All dependence on the angular coordinates is included in spherical harmonics 
which form an irreducible representation of the group 0(3), and the dependence on the 
single remaining radial coordinate is governed by a one-dimensional Schrodinger equation. 
Similarly, the conformal expansion of distribution amplitudes in QCD aims to separate 
longitudinal degrees of freedom from transverse ones. All dependence on the longitudinal 
momentum fractions is described by orthogonal polynomials that form an irreducible rep- 
resentation of the so-called collinear subgroup of the conformal group, SL(2,M), describing 
Mobius transformations on the light-cone. The transverse-momentum dependence (scale- 
dependence) is governed by simple renormalisation group equations: the different partial 
waves, labelled by different "conformal spins", behave independently and do not mix with 
each other. Since the conformal invariance of QCD is broken by quantum corrections, 
mixing between different terms of the conformal expansion is absent only to leading log- 
arithmic accuracy. Still, conformal spin is a good quantum number in hard processes, up 
to small corrections of order a^. The application of conformal symmetry to the study 
of exclusive processes to leading twist has become a vast field whose current status is 
reviewed in Ref. [11]. 

Since quark mass terms break the conformal symmetry of the QCD Lagrangian explic- 
itly one might, at first glance, expect difficulties to incorporate SU(3) breaking corrections 
into the formalism. In fact, however, the inclusion of quark mass corrections is straight- 
forward and generates two types of effects. First, matrix elements of conformal operators 
are modified and in general do not have the symmetry of the massless theory. This is 
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not a "problem", since the conformal expansion is designed to simplify the transverse 
momentum dependence of the wave functions by relating it to the scale dependence of the 
relevant operators. This dependence is given by operator anomalous dimensions which are 
not affected by quark masses, provided they are smaller than the scales involved. Second, 
new higher twist operators arise, in which quark masses multiply operators of lower twist. 
These additional operators generate contributions to higher twist DAs, which have been 
discussed in [9] and are irrelevant for the present investigation. 

The conformal expansion of DAs is especially simple when each constituent field has 
a fixed (Lorentz) spin projection onto the light-cone. In this case, its conformal spin is 

j^lil + s), (2.3) 

where / is the canonical dimension and s the (Lorentz) spin projection. Multi-particle 
states built of constituent fields can be expanded in increasing conformal spin: the lowest 
possible spin equals the sum of spins of the constituents, and its "wave function" is given 
by the product of one-particle states. This state is nondegenerate and cannot mix with 
other states because of conformal symmetry. Its evolution is given by a simple renormali- 
sation group equation and the corresponding anomalous dimension is the smallest one in 
the spectrum. This state is the only one to survive in the limit — > oo and is usually 
referred to as "asymptotic distribution amplitude" . 

As for the leading- twist quark- antiquark DAs we are interested in, the partial wave 
expansion reads 

oo 

(j){u) = 6uu anC^/2(2M - 1), (2.4) 

n=0 

where Cn'^{2u — 1) are Gegenbauer polynomials. The dimension of quark fields is Z = 3/2 
and the leading twist distribution corresponds to positive spin projection s — +1/2 for 
both the quark and the antiquark. Thus, according to (2.3), the conformal spin of each 
field is jq = jg = 1; the asymptotic distribution amplitude is 6uu. The Gegenbauer poly- 
nomials correspond to contributions with higher conformal spin j + n and are orthogonal 
over the weight function 6uu. 

Note that ao — 1 due to the normahsation condition (2.2). In the hmit of exactly 
massless quarks only terms with even n survive in Eq. (2.4) because of G-parity invari- 
ance. The conformal expansion, however, can be performed at the operator level and is 
disconnected from particular symmetries of states such as G-parity. The expansion (2.4) 
is, therefore, valid for arbitrary n, and it is precisely the odd contributions to the expan- 
sion which, being proportional to manifestly SU(3) breaking effects like the difference of 
quark masses, induce the most tangible SU(3) breaking effects. 

As mentioned before, conformal invariance implies that partial waves with different 
conformal spin do not mix under renormalisation to leading-order accuracy, which means 
that the Gegenbauer moments a„ in (2.4) renormalise multiplicatively: 

«.(/.) =a.(/.o)^^J (2.5) 
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with (3q — 11 — (2/3)n/. The one-loop anomalous dimensions are [12] 



/ 2 n+l \ 

= ^(") = ^-['- (n+l)(n + 2) + ' g '/'J ' 

7fn) = C'^(^l + 4^1/jj . (2.6) 

7(0) is the anomalous dimension of the local current and vanishes for vector and axialvector 
currents. Note that the DAs of the pseudoscalar K and the longitudinally polarised K* 
have the same anomalous dimensions. As the 7's are positive and increase with n, the 
effect of running to large scales is to damp the effects of higher-order Gegenbauers, so 
that the DAs approach their asymptotic shape. Prom the above-stated scaling-properties, 
it is evident that a discussion of the shape of DAs is most conveniently done in terms of 
Gegenbauer moments. 



3 QCD Sum Rules for Moments 
3.1 QCD Sum Rules 

As QCD sum rules constitute an estabhshed method for the calculation of static hadro- 
nic matrix elements like decay constants (as opposed to dynamical quantities like form 
factors), we refrain from delving into a comprehensive explanation of technicalities, for 
which we refer to the original papers [13] and to recent reviews [14] instead. The key 
feature of the method is the use of analyticity to relate the local short-distance operator 
product expansion (OPE) of a correlation function of two currents, 

n = i /" d'ye"^y{0\TMy)MO)\0) = ^ (:7„(g2)(0„) = n°™ (3.1) 

around y = (as opposed to a light-cone expansion aound = 0, which is appropriate 
for form factor calculations, cf. [10]) valid for = —q^ -C 0, to its dispersion relation in 
terms of hadronic contributions. 



poo 

^0 



ds P^'' ^ n'^^'i, (3.2) 



where p{s) is the spectral density of the correlation function along its physical cut. The 
OPE yields a series of local operators of increasing dimension whose expectation values 
iPn) in the nonperturbative (physical) vacuum are the so-called condensates. In the sum 
rules analysed in this paper, we take into account the condensates listed in Tab. 1. The 
representation of the correlation function in terms of hadronic matrix elements can be 
written as 
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(qq) = (-0.24 ± 0.01)3 GeV^ 
G^) = (0.012 ± 0.006) GeV^ 
(qagGq) = (0.8 ± 0.1) GeV^ {qq) 


{ss) = {0.8±0.1){qq) 
{sagGs) = (0.8 ± 0.1){qagGq) 


m,(2GeV) = (110±20)MeV < — > m,(l GeV) = (149 ± 27) MeV 
a,(lGeV) - 0.513 i — > Aqct ° " 372 MeV 



Table 1: Input parameters for sum rules at the renormalisation scale = 1 GcV. The 
value of nis is obtained from quenched lattice calculations as summarised in [15]. 

where uim is the mass of the lowest-lying state couphng to the currents Ji,2 and p'^°^^ 
parametrises all contributions to the correlation function apart from the ground state. /, 
the residue of the ground state pole is the quantity one wants to determine. A QCD sum 
rule that allows one to do so is obtained by equating the representations (3.1) and (3.2) 
and implementing the following (model) assumptions: 

• p^°'^* is approximated by the spectral density obtained from the OPE above a certain 
threshold, i.e. p^"^^ p'^^^{s)9{s — sq) with sq ~ {mM + A)^ being the continuum 
threshold, where A ~ O(Aqcd) is an excitation energy to be determined within the 
method. This assumption relies on the validity of global quark-hadron duality; 

• instead of the weight-functions l/(g^)" and l/{s — g^), one uses different weight- 
functions which are optimised to (exponentially) suppress effects of p{s) for large 
values of s and at the same time also suppress high-dimensional condensates by 
factorials. This is achieved by Borel transforming the correlation function: Bl/{s — 
q^) = 1/M^ exp(— s/M^). A window of viable values of the Borel parameter 
and the continuum threshold Sq has to be determined within the method itself by 
looking for a maximum region of minimum sensitivity (a plateau) in both and 
■So; 

• the OPE of n can be truncated after a few terms. As we shall see (and as is well 
known), this condition is fulfilled only for low Gegenbauer moments. 

After subtraction of the integral over p^^^ above Sq from both sides, the final sum rule 

reads 

B^^^U^^^ ^±-^£dse-^"^'p^^^{s) = ^e-M/M^, (3.3) 

which gives the hadronic quantity / as a function of the Borel parameter and the 
continuum threshold Sq (and the condensates and short-distance parameters from the 
OPE). 

Even Gegenbauer moments can be determined from diagonal correlation functions of 
type 

I [ d'ye^^y{0\Tq{y)Ts{y)mmMm, (3-4) 
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with a suitably chosen Dirac structure F. This form of the correlation function differs 
from Refs. [5, 7], where local operators with an arbitrary number of covariant derivatives 
were used rather than the nonlocal operator s(0)r[0, z]q{z). We find the calculation with 
nonlocal operators very convenient, as it allows one to calculate all moments in one go. 
Specifying for instance to ^1*, the sum rule reads 

BsubHOP^ = (/i)'e-^*/^^ ^^'di.e^^(''-)4(«), (3.5) 

where also n*-*^^ is expressed as integral over u, which naturally emerges as Feynman 
parameter in the calculation, and comes with the same weight function exp[iu[q ■ z)). 
Sum rules for individual Gegenbaucr moments are obtained by expanding both sides in 
powers of {q ■ z), or effectively replacing 

^iui,.) ^ cr{2u - 1), j\ue^-^^-^i{u) ^ '^^X{2n^t)'^^ 

For odd moments, one analyses nondiagonal correlation functions of type 

i j dV^^(0|r^(y)F2s(y)s(0)Fi[0,z]g(z)|0) (3.6) 

with structures Fi and F2 of opposite chirality [8, 5]. The Fj appropriate for K*^ and 
K]_ are given in App. A, where we also give complete expressions for perturbative and 
quark-condensate contributions to 0(q;s), as well as tree-expressions for the dimension 5 
mixed condensate. 

Formulas for the correlation functions and Gegenbauer moments are collected in the 
appendices. Note in particular the scaling of different contributions to the Gegenbauer 
moments in n, the order of the moment: nonperturbative terms increase with positive 
powers of n with respect to the perturbative contribution. For large n, this behaviour 
upsets the usual hierarchy of contributions to the OPE, where nonperturbative terms 
are expected to be a moderately sized correction to the leading term. The origin of this 
behaviour can be easily understood from the fact that in the local expansion the vacuum 
fields have exactly zero momentum, which yields the ^-function terms in (A. 3) and (A. 4). 
This amounts to a multipole expansion of the DA around its endpoints, which is justified 
if one is only interested in gross features of the DA, like the first few moments, but which 
is clearly inappropriate for extracting more detailed information on the shape. In the 
present paper we adopt the viewpoint that at least the first two Gegenbauer moments 
can be reliably obtained from local sum rules of the type discussed above. 

3.2 A Short Discussion of Decay Constants and Lattice Results 

The decay constant of the pseudoscalar K meson is well known from — > /x"'"i/^(-|-7) 
and quoted as [16] 

fK = (159.8 ± 1.3) MeV. 
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The decay constant can be extracted from the branching ratio of r — > K* via 




With \ Vus\ = 0.220 and the other parameters taken from [16], one finds 

fl = (217 ± 5) MeV. 



It is interesting to compare this value with a "postdiction" from QCD sum rules. We use 
the sum rule quoted as (C.4) in [9] with the input parameters as in Tab. 1. Prom the 
result plotted in Fig. 1 we conclude 



and So = (1.8 ± 0.1) GeV^ The plateau in extends from 1 to 2 GeV^. Although 
the error only includes the impact of varying the input parameters of Tab. 1 within their 
respective ranges and does not account for any systematic uncertainties, the agreement 
with the experimental result is remarkably good. 

The value of the third constant, f^, has not been measured yet, but has been calculated 
in quenched lattice QCD and from QCD sum rules. The lattice result for the ratio of the 
K* decay constants is [17] 



using the formulas (C.4) and (C.5) in [9], which include NLO radiative corrections. The 
optimum continuum threshold for turns out to be sq = (1.2 ± 0.1) GeV^. The fact 
that So for is smaller than for K*^ is in agreement with the findings of Ref. [7] and 
related to the fact that, with the Dirac structure T chosen as specified in App. B, the 
sum rule does not only receive contributions from K*, but also from the ground state in 
the opposite parity channel, the l"*" state i^i(1270). This contribution has to be included 
in the continuum, which results in a low value of sq. Note that this is not an external 
condition imposed by us when evaluating the sum rule, but emerges naturally from the 
criterion of stability and the presence of a plateau in M^. We consider the ratio of 
decay constants rather than the constants themselves, as effects of unknown higher order 
corrections are expected to cancel in the ratio. The sum rules are evaluated at the low 
scale 1 GeV and the scaling from 1 to 2 GeV is done using NLO renormalisation group 
improvement. Note that the dependence on the Borel parameter largely cancels in the 
ratio of decay constants and also the dependence on sq is rather mild. The agreement 
with the lattice results is remarkably good. Using the experimental value for /|-, QCD 
sum rules hence predict, to NLO accuracy and including 2-loop running. 



fHsR) _ ^225 ± 7) MeV 



Prom QCD sum rules, on the other hand, we find, cf. Fig. 2, 






0.74(03), 



/^(l GeV) = (182 ± 10) MeV, 

/i^(4.8 GeV) = (149 ± 8) MeV. 
This updates the result obtained in Ref. [9] 



/^(2 GeV) = (160 ± 9) MeV, 



(3.7) 
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0.205 I — ' ' ' ' ' ' ' ^ 

0.6 0.8 1 1.2 1.4 1.6 1.8 2 

M^[GeV]^ 

Figure 1: ff^ from QCD sum rules as function of the Borel parameter. The spread between 
the two curves corresponds to the uncertainty induced by varying the input parameters 
within their error margins. 
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o 




0.6 0.8 



1.2 1.4 1.6 1.8 



M^[GeV]^ 



Figure 2: The ratio /^(l GeV)//|- as function of the Borel parameter M^. Interpretation 
of curves as in previous figure. Sq is varied between 1.1 and 1.3 GeV^ for and 1.7 and 
1.9 GeV^ for /|. 



3.3 The Gegenbauer Moment ai 

Let us now calculate the first moments of the K and K* DAs. Previous determinations 
go back to Refs. [5, 9], where the following values were obtained: 

ai(lGeV) = 0.17, af(l GeV) = 0.19 ± 0.05, a^(l GeV) = 0.20 ± 0.05. 

These results are valid at LO in QCD and were obtained using QCD sum rules which, 
unfortunately, partly suffered from mistakes. In particular, we find a different sign of 
the perturbative contribution w.r.t. the formulas given in Ref. [5, 9]. We have carefully 
checked that the sign we obtain is indeed the correct one: for the nth Gegenbauer moment, 
the leading order perturbative contribution is the same for all three correlation functions 
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-0.05 I < < < < < < < r 

-0.1 

-0.15 - 




-0.45 I — ' ' ' ' ' ' ' ^ 

0.6 0.8 1 1.2 1.4 1.6 1.8 2 

M^[GeV]^ 

Figure 3: ai(lGeV) as function of the Borel parameter M^. The spread between the two 
curves corresponds to the uncertainty induced by varying the input parameters within 
their error margins. Sq is fixed at 1.8 GeV^. 

-0.15 r-^ < < < < < < ^ 




0.6 0.8 1 1.2 1.4 1.6 1.8 2 
M^[GeV]^ 

Figure 4: af(lGeV); notations and conventions as above. 



-0.1 r-^ < < < < < < r 




0.6 0.8 1 1.2 1.4 1.6 1.8 2 
M^[GeV]^ 



Figure 5: aj^{lGeV); So is fixed at 1.2 GeV^; other notations and conventions as above. 
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and is given by 

It appears that the factor (—1)" has been missed in [5]. Due to the change in sign of 
the perturbative contribution, the leading order sum rule becomes numerically unstable: 
perturbation theory and the mixed condensate yield negative contributions, the quark 
condensate a positive one, and the sum is close to zero. This is in stark contrast to the 
statement made in Ref. [5], according to which the sum rule is dominated by the quark 
condensate contribution and yields a positive result for all ai. 

The different signs of individual contributions, together with the fact that none of them 
is numerically dominant, entails that no meaningful result can be extracted from the lead- 
ing order sum rule. The situation improves, however, if one includes radiative corrections. 
We have calculated 0{as) corrections to both the perturbative and the quark condensate 
contribution and find that they reduce the size of the quark condensate contribution, but 
increase that of the perturbative contribution. As a result, the perturbative contribution 
becomes numerically dominant, as is actually expected for a "good" sum rule. By varying 
all the input parameters of Tab. 1 within their respective ranges we obtain 

ai(lGeV) = -0.18±0.09, a!(l GeV) = -0.4±0.2, a^(l GeV) = -0.34±0.18. (3.8) 

The Borel window is taken to be 1 to 2 GeV^, as motivated by the results of the previous 
subsection. Note that the l.h.s. of the sum rules for af'"*" contains the factor which 
we substitute by their respective sum rules instead of using the values determined in the 
previous subsection, the reason being that one expects unknown higher order corrections 
to cancel in the ratio. The dependence of the individual sum rules on the input parameters 
is plotted in Figs. 3, 4 and 5. The error bars in (3.8) are rather large, which is due to 
the fact that the first Gegenbauer moments are explicitly proportional to SU(3) breaking 
quantities, i.e. or the difference of condensates, e.g. {qq) — {ss), which come with a 
considerable uncertainty. 

We would also like to add that the use of nondiagonal sum rules has been met with 
criticism. It has been argued in Ref. [18] that these sum rules may suffer from large 
contributions of higher resonances or, in the case of the K, from instanton contributions to 
the pseudoscalar current. We can actually estimate the amount of possible contamination 
of these sum rules by studying their local limit, yielding ao, which is 1 by definition. 
In Fig. 6 we plot the sum rules for oq for both the K and the K* (there is only one 
sum rule for the K*), for central values of the input parameters, Sq — 1.2 GeV^ for K*, 
Sq — 1.8 GeV^ for the K, and replacing /^'" in the denominator by their respective sum 
rules. The plot shows no sizeable contamination for the K, and a moderate one for the 
K*; both results lie in the ball-park of the expected accuracy of QCD sum rules and 
come with uncertainties from the input parameters, in particular for ao of the K, which is 
directly proportional to m^. Radiative corrections are important in both cases and drag 
the result closer to ao = 1. This result strengthens our confidence in the suitabihty of 
nondiagonal sum rules for extracting meaningful values of ai. 



11 



1.3 
1.25 

1.2 
1.15 

1.1 
1.05 
1 

1 1.2 1.4 1.6 1.8 2 
M^[GeV]^ 

Figure 6: gq from nondiagonal sum rules for K (lower curve) and K* (upper curve) for 
central values of input parameters. 

3.4 The Gegenbauer Moment a2 

Sum rules for even moments and the effects of SU(3) breaking have been studied in var- 
ious papers, e.g. [8, 5, 9]. None of these papers, however, does contain a complete set of 
sum rules for Gegenbauer moments. Moreover, while checking our two-loop and nonper- 
turbative calculations against the formulas collected in Ref. [9], we found that the terms 
in irisisagGs), given correctly in Ref. [5] for the moments of K^*\ have not been correctly 
translated into Gegenbauer moments in Ref. [9]; we have calculated these contributions 
anew for all three DAs, using nonlocal currents, and confirm the results of Ref. [5]. We 
have also recalculated the four-quark contributions and confirm the expressions quoted in 
[9]. The complete set of (hopefully) correct sum rules is given in App. B. 
With the input parameters from Tab. 1 we obtain 

02(1 GeV) = 0.16 ±0.10, 4(1 GeV) = 0.09 ±0.05, a^(l GeV) = 0.13 ± 0.08 (3.9) 

by varying all the input parameters of Tab. 1 within their respective ranges. These 
numbers have to be compared with those quoted in the first reference in [10] and in 
Ref. [9]: 

02(1 GeV) = 0.2, 4(1 GeV) = 0.06 ±0.06, a^(l GeV) = 0.04 ± 0.04. 

The discrepancy to (3.9) is mainly due to the wrong expressions for the contribution of 
the (sagGs) condensate used in [9]. 

Like for ai, the relative errors of the moments (3.9) are considerably larger than the 
ones quoted for the decay constants. The reason is the absence of a stability plateau 
in as the perturbative contribution is of 0{as) and small. For the extraction of the 
Gegenbauer moments we thus have to rely on the optimum values of M^, 1 to 2GeV^, 
and So determined from the sum rules for the decay constants. 

3.5 Models for Distribution Amplitudes 

Nearly all models for DAs rely on a truncated conformal expansion, which is not too 
surprising in view of the fact that it is only moments (or Gegenbauer moments) that are 
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0.05 ^ ' ' ' ' ' ^ 

0.8 1 1.2 1.4 1.6 1.8 2 
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Figure 7: a2(lGeV) as function of the Borel parameter M^. The spread between the two 
curves corresponds to the uncertainty induced by varying the input parameters within 
their error margins. Sq is fixed at 1.8 GeV^. 
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Figure 8: a2(lGeV); notations and conventions as above. 
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Figure 9: a^(lGeV); Sq is fixed at 1.2 GeV^; other notations and conventions as above. 
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K Kl Kl 


ai(lGeV) 
ai(2GeV) 
ai(4.8GeV) 


-0.18 ±0.09 -0.40 ±0.20 -0.34 ±0.18 
-0.15 ±0.07 -0.32 ±0.16 -0.28 ±0.14 
-0.13 ±0.06 -0.28 ±0.14 -0.23 ±0.11 


02(1 GeV) 
02(2 GeV) 
02(4.8 GeV) 


0.16 ±0.10 0.09 ±0.05 0.13 ±0.08 
0.11 ±0.07 0.06 ±0.04 0.10 ±0.06 
0.09 ±0.06 0.05 ±0.03 0.08 ±0.05 



Table 2: First and second Gegenbauer moments ai and 02 of the leading-twist DAs of K 
and K* mesons evaluated at different scales. 



accessible to direct calculation to date. Truncation after the first few moments implies 
that only gross features of the DAs are taken into account, which is sufficient if physical 
amplitudes are obtained by convoluting with a smooth perturbative scattering amplitude 
- or if the characteristic scale of the process is large enough for the logarithmic damping 
of high-order Gegenbauer moments to have become effective. In this paper we assume 
that this is indeed the case and that a truncation after n = 2 yields a meaningful approx- 
imation. We then obtain the following model for the three leading-twist DAs of strange 
mesons containing an s quark: 



(f)K{u, IJ>) — 6u{l — u) 



1 + 3ai(ii)(2u - 1) + a2(/i) - (5(2u - if - 1) 



(3.10) 



To obtain the DAs for mesons with an s antiquark, one has to replace u \ — u. The 
values of the Gegenbauer moments a\ and 02 are calculated from QCD sum rules at the 
low scale 1 GeV and scaled up to 2 and 4.8 GeV, to LO accuracy;^ numerical values are 
collected in Tab. 2. Eq. (3.10) is ready for use in phenomenological applications. 



4 Summary and Conclusions 

We have presented a comprehensive study of SU(3) breaking corrections to the decay 
constants and the leading-twist light-cone distribution amplitudes of K and K* mesons 
from QCD sum rules. We have corrected a few computational mistakes in previous works 
and included radiative corrections to the leading contributions to odd moments of the 
distribution amplitudes. Our main results are summarised in Eqs. (3.7), (3.10) and Tab. 2. 
We find in particular that the first Gegenbauer moments arc negative for all three DAs, 
in contrast to previous determinations. The decay constants ]k and ]\ can be extracted 
from experiment with good accuracy, whereas has to be determined theoretically. 
Its accuracy is limited mainly by the uncertainty of the input parameters, whereas the 
accuracy of the Gegenbauer moment a\ is in addition limited by the lack of a proper 
plateau region in tlie Borel ])aiaincter, which suggests that the uncertainty is dominated 

^NLO scaling mixes all Gegenbauer moments, so that to this accuracy the model (3.10) receives 
radiatively generated corrections in a„>2. 
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by systematics. The latter limitation also affects 02- This implies that further refinement 
of the sum rules by calculating even higher order corrections is not likely to improve their 
accuracy. Any decisive improvement has to come from a different method, for which 
lattice QCD appears to be a natural candidate. It is to be hoped that the preliminary 
results for the second moment of the twist-2 pion DA reported in [4] will be improved and 
also extended to the K and other mesons. 

Eq. (3.10) and Tab. 2 present an approximation to the full DA to NLO in the conformal 
expansion. A question we have not touched upon in this paper is the suitability of this 
approximation for calculating physical amplitudes of type 



where T{u) is a perturbativc scattering amplitude. Obviously the answer to this question 
depends on the functional dependence of T on u. The standard argument to justify the 
validity of a truncated conformal expansion is that for sufficiently smooth functions T{u) 
higher order Gegenbauer polynomials in ^{u) with their highly oscillatory behaviour are 
effectively washed out. It is also argued that, for large scales, higher order Gegenbauer 
moments are suppressed by renormalisation group scaling, so that for most applications 
the asymptotic DA is sufficient. These arguments are usually presented in a qualitative 
rather than a quantitative way and it would be interesting to study, within a well-defined 
model that does not rely on conformal expansion, how sensitive physical amplitudes ac- 
tually are to a truncation of the conformal expansion, and at what scales logarithmic 
damping becomes effective. These questions will be addressed in a future publication. 
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Appendix 

A Results for Nondiagonal Correlation Functions 



where ri^2 are Lorentz-structures suitable for the projection onto the meson M, M G 
{K, K]_, Ku}, = is a light-hke vector and 



the path-ordered gauge-link joining the quark fields in the bilinear si^^Yiqi^z). The above 
sign-convention for g implies that the covariant derivative is given by = — igA^. 



■phys 




We calculate the correlation function 




(A.l) 
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r2 


projection V : IIm VUm 


K 




75 


-1/iq-z) 




z 




i/{q-zf 


Kl 




la 





Table A: Dirac structures used in Eq. (A.l). The choice of ri^2 together with the projection 
V projects onto the twist-2 structure and ehminates factors i{q ■ z) picked up by the 
traces. D is the number of dimensions in dimensional regularisation. We use the notation 
a — a^'y^ for an arbitrary 4- vector a. 



Odd Gegcnbaucr moments are most conveniently determined from nondiagonal correlation 
functions, involving one chiral-odd and one chiral-even current. To be specific, we choose 
the currents given in Tab. A and perform an operator product expansion of IIm up to 
dimension 5 condensates. We express the perturbative results in the form 



1 

■v2 



p-QPcit.th ^ — rns{ij) In — |- j e*"*^^'^-* u -KMiu; /j.) + terms analytic in g^, 



where u is the momentum fraction carried by the s quark in the K or K* meson; u = 1—u 
is the momentum fraction carried by the u or d antiquark. Note that one has to define 
u in such a way that the exponential reads exp{iu{q ■ z)), rather than exp{iu{q ■ z)), as 
it has to match the corresponding factor in the hadronic parametrisation of the corre- 
lation function, cf. Eq. (3.5). The 0{as) corrections are new, and we find the following 
expressions: 

^ ^ a, (51-2n^ , -q^ / 9 , \ 2 (-3 + 4:u) Inu In^ii 
nK = l + CF-^< + ln^ ---Inw + ^ ^ 



47r|^3 fji^ \ 2 J u u 

+ (-3-2 Inn) lnu + 2ln^u'^ + 4L2(u) + 2L2(--) -2- (L2(u)-L2( — )) }, 

a f 33 — 27r^ /5 \ —a^ 

TT;^* = 1 + i V2\n^u- {'i + 2\uu)\uu- - H-lnw In ^ 



47r [ 3 ' ' \2 ) [i^ 

^ ' hiu = \ ' ^ +4L2(^) + 2L2(-m/«) \ , 



u u u u 



a,r33 - 27r2 2- o2-«, l + w,2 (o^^'^A^ ''f 
ttk* — 1 + Cp — s h 2 m u — 6 — z — z — In u — I 6-\ — — In 



'-^ 47r [ 3 u u \ u J jj? 

- (^3 + - 1^ ^2(ii) + ^L2{u) + 2L2{-u/u) I . (A.2) 
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r2 


projection V : IIm VUm 


K 






l/{q.zf 


Kl 


z 


z 


l/{q-zf 


Kl 




ai'Pzp 


1 1 

2 - D (g ■ 



Table B: Dirac structures and projections used for calculating even moments from diagonal 
correlation functions 



A collection of loop integrals necessary to perform these calculations can be found in 
App. C. We have checked that the above expressions reproduce the correct anomalous 
dimensions for the Gegenbauer moments, 7''^-^ and l[n)-> Eq- (2.6). We have also checked 
that in the local limit, i.e. ^ — >• 0, e^\){iu{q- z)) — >• 1, the expressions for K^^ and K]_ agree: 
VH-K* = VH-K*^. Note that the sign of the leading term differs with respect to [8, 5]. 
We have also calculated 0{as) corrections to the condensate terms and find 



Q Jo 



{ss){Siu) + ^Cpi^6iu) (^51n-|- 



+ (1 - 5m,\\) (2 - 36m,± + {3Sm,± - 4) In -^^ ^-2u{l- 5m,j 



X 1 + 25m,\\ + ln{uu) + In 



+ 2 



M / , / _x , —q 
— 2 — ln(uu) — m — - 

u \ fJ,^ 



+{qq)(u ^ u) 
where the [ ]+ prescription is defined as 



[fiu)U^fiu)-5iuo) f'dvfiv), 

Jo 



(A.3) 



if / has a simple pole (modulo logarithms) at < Uq < 1. 

For the contribution from the dimension 5 mixed condensates (qagGq) and (sagGs), 
we restrict ourselves to the tree-level contribution and find 

^^{mog) du^^ii-^^ [{sagGs){{5M,\\ + 5m,a.)5{u) + 5\u)} + {qagGq){u ^ u}] . 

•J? Jo 

(A.4) 



B Sum Rules for Moments 

The Dirac structures and projections used for calculating the diagonal correlation function 
Eq. (3.4) are collected in Tab. B. 
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The nonperturbative corrections to the sum rule for even moments of the K can be 
extracted from [5]. We have recalculated the SU(3) breaking terms explicitly and find 
agreement with Rcf. [5]. The perturbativc terms can be extracted from Ref. [19]. We have 
also recalculated these terms in the nonlocal operator formalism and confirm the result 
quoted in [19]. The complete sum rule for even Gegenbauer moments of the K reads 

2(2n + 3) ■'^ "^^^ 27r2 TT V / 

duuuC'J\2u - 1) In^ ^ + (n + l)(n + 2) 



X 



+ ^ ^ (» + 1)(„ + 2) + 1^ ^^^^^ (» + mn + 2)^ (B.1) 

We have also checked the nonperturbative terms in [9] and could not confirm the terms 
in TTisisagGs) . We thus find it appropriate to present here the (hopefully) correct sum 
rules: 



27r2 



^ (l - e-'o/M') duuuCT{2u - 1) W 



+ ^ ^s{ss){n + l)(n + 2) + ^ G^) (n + l)(n + 2) 

^ m,(sa5G's)n(n + l)(n + 2)(n + 3) - ^^^^^ (gg)(ss)(n + l)(n + 2) 



((5^5)' + (^"^)')(^ + + 2)^ (B.2) 



2(2n + 3) 

1 a, 
2^ 



^ (^1 - e-^oVM^ j j\uuuC^J\2u - 1) (inw + Inw + In^ ^) 



+ 24^ G^) (n^ + 3n - 2) - ^ m,(sayGs)(n + l)(n + 2){n' + 3n + 4) 
+ 8171^ + ^^^^'^^"^ " ^^^"^ + ^^^"^ + ^^^"^ + + 2M^ ""^^^^^^^ 
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Sum rules for odd moments are obtained from the formulas given in App. A by Borel 
transforming and replacing exp(m(g • z)) Cn'^{2u — 1). In this way, and using the 
properties of Gegenbauer polynomials as for instance collected in [20], which amount to 
the following replacements: 

5{u) ^ (-1)" 1 (n + l)(n + 2), 5{u) ^ i (n + l)(n + 2), 

5\u) (-1)" ^ n{n + l){n + 2){n + 3), S'{u) ^^n{n + l){n + 2)(n + 3), 

we find for the 0(q;°) contributions for odd n: 

BsubPni,") = -A-^ms{l- e-°/^^) + ^ (n + l)(n + 2)((.-.) - (qq)) 

+ 12^ + + 2)(n + 3){{qagGq) - {sagGs)) , 
SsubPnJ] = - A (^1 _ e--/*^^) + _1_ + + 2){{ss) - {qq)) 

+ in + iTin + 2)\{qagGq) - {sagGs)) , 

-BsubPnS^I = B^^^VIi^ . (B.4) 

Note that the sign for the leading order perturbative contribution is different with re- 
spect to (C.9) and (C.IO) in Ref. [9] and also with respect to [8, 5]. We have carefully 
checked that we obtain the correct sign in the local limit z ^ Q and that our result 
stays unchanged for a different choice of the position of the nonlocal current in coordinate 
space, for instance s(0)ri[0, z\q{z) — >• s{—z)Vi[—z, z\q{z). As it turns out, also the mixed 
condensate-contributions differ from the formulas given in [9] . 

To obtain the 0(as) corrected sum rules for odd moments from the formulas given in 
App. A, we also need the following continuum-subtracted Borel transforms: 



-Bsub In 











-e- 


iO 






-e- 


iO 





= - (l - , 



1 , -(/^ - iO If , „. / So \ 1 



which complete the set of formulas needed to translate the correlation functions obtained 
in the previous appendix into QCD sum rules for odd moments. 
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Figure 10: 0{gg) and 0{gl) Feynman rules induced by the nonlocal current in (A.l). 

C Loop Integrals 

To the benefit of apprentices, and also for future reference, we collect in this appendix 
relevant one- and two-loop integrals. The main difference as compared to usual calcula- 
tions is the nonlocal vertex induced by s{0)Ti[0, z]q{z), which, to order a^, gives rise to 
the Feynman rules shown in Fig. 10. 

The master one- loop integral is given by {z^ — 0, D — A + 2e) : 



-i(k—q)z 



{kzf 



X f e^«(«-^) (C.l) 
Jo 



where [rf^fc] i/(47r)^/2 ^ d^k/{2n)^. 

As for two-loop integrals, one needs for instance 



kH'^{k-iy{k-qf r(l-e) 

X— ^ f c^iie^"(^-") / dyy'y-^+' ((1 - uy^-^ - u 

Al ' z) Jo t Jo 



a-l> 



r(l + 2e) r(6 + l + 2e) / 

Generally, all loop-integrals with additional exponentials can be calculated conveniently 
using Feynman parameters. The calculation is further simplified by the fact that one only 
needs the imaginary part in g^, which implies that finite integrals need not be calculated. 
"Overlapping exponentials" like for instance 



dx 
Jo Jo 



dy ex.p{i{x + yx){q ■ z)) f{x,y) 

Jo Jo 

can be rewritten as 



/ du exp(iu{q ■ z)) dx dy 5{x + yx — u)f{x,y), 
Jo Jo Jo 



which allows one to reduce all contributions to the canonical form. 
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